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1. Introduction 

In addition to being an important tool for modeling nucleus-nucleus collisions (AA), proton- 
nucleus (pA) collisions are also interesting by themselves as they probe nuclear effects related to 
multiple scattering and possible hadronization in nuclear matter. The origin of these effects is still 
under debate. In recent years most of theoretical activity has focused on novel high-energy nuclear 
effects, such as partem saturation, yet it seems that past and present energies (RHIC) do not provide 
enough phase space for these effects to appear. It will be our task in this paper to review some 
well-known low-energy models and to introduce high-energy effects (shadowing) to describe data 
on light and heavy particles in pA collisions in the energy range y/s = 17.3 — 5500 GeV. 

Our starting point is noticing that a significant change in the underlying dynamics of a hadron- 
nucleus collision takes place with growing energy of the incoming particles. At low energies, the 
total cross section is well described within the probabilistic Glauber model |jlj|. In the reggeon 
approach these interactions are described by so-called planar diagrams depicted in Fig. [I] (left). At 
higher energies, E > Eq ~ m^fLRA (jU is a characteristic hadronic scale, fL ~ 1 GeV, and Ra is the 
radius of the nucleus) corresponding to a coherence length 

1 

(1.1) 



2m^x 



the typical hadronic fluctuation length can become of the order of, or even bigger than, the nuclear 
radius and there will be coherent interaction of hadron constituents with several nucleons of the 
nucleus. In this energy range, the contribution from planar diagrams is damped by a factor ~ l/E 
[Q] and the dominant contribution arises from non-planar diagrams, shown in Fig. [T] (right). The 
sum of all diagrams was calculated by Gribov ^ |J], who introduced a correction to the Glauber 
series by taking into account the diffractive intermediate states. The forward hadron-nucleus QiA 
or y*A) scattering amplitude can then be written as the sum of diagrams shown in Fig. ||, i.e. as 

o f A = Aa rN + afil + ... , (1.2) 



In Eq. (1.2), the first term simply represents the Glauber elastic contribution and subsequent terms 
describe multiple interactions of the incoming probe with the nucleons in the target nucleus. The 
space-time picture analogy to the Glauber series is however lost, as the interactions with different 
nucleons of the nucleus occur nearly simultaneous in time. This phenomenon is related to inelastic 
shadowing corrections. 

An additional effect which comes into play at high energies, is the possibility of interactions 
between soft partons of the different nucleons in the nucleus. In the Glauber-Gribov theory this 



Figure 1: Planar (left) and non-planar (right) diagram of multiple scattering. The former decreases as ~ l/E 
above the critical energy, the latter is controlled by f„„„ effects. 
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corresponds to interactions between Pomerons. The necessity to include such "enhanced" diagrams 
at high energies can be related to unitarization of all amplitudes. In particular, diagrams involving 
triple-Pomeron interactions are related to large-mass intermediate states in Fig. |2] which give a 
dominant contribution to shadowing. There is a connection between these effects and saturation 
effects in the parton picture. 

In pA collisions, nuclear effects are usually discussed in terms of the power-law parameteri- 
zation 

dp dp 

where o p a (<J p n) i s the inclusive cross section off a nucleus (nucleon). The function a(xp ) char- 
acterizes nuclear effects at different longitudinal momentum fractions of the produced particle. 
For a large range of energies, a exhibits a very interesting scaling in xp for both light and heavy 
particle production as seen in Fig. ||. For J/yr, <x(xp) decreases from 0.95 at xp « to values 
a ~ 0.75 at xp ~ 0.8 thus indicating an increase of absorption as xp increases No scaling in 

Bjorken x of the nucleus, X2, is observed, indicating breaking of QCD factorization at these ener- 
gies. The suppression of light particles follow a similar trend. Another striking feature is the rather 
large suppression of low-energy data at xp = 0. These features are reproduced by models invoking 
mechanisms of attenuation and/or energy-loss in nuclear matter [p|-[TT|]. 

Recent data from RHIC experiment on charged hadron [|12|] and J/\jf [13] production in dAu 



collisions at y^s = 200 GeV are also shown in Fig. || (data from [12] have been integrated over p± 
between 0.7-1.0 GeV/c). The suppression at xp = is smaller compared to lower energies contrary 
to the expectation of many theoretical models. Despite the limited kinematics accessible in collider 
experiments, the RHIC data also hints towards a breaking of xp scaling. 

The behavior of a(xp) with energy allows a natural explanation within the Gribov theory of 
multiparticle production described above. At very high energies Abramovsky-Gribov-Kancheli 



(AGK) cutting rules [ |15| ] lead to cancellation of the Glauber-type diagrams in the central rapidity 
region and only enhanced diagrams, or in other words shadowing of small-x partons, contribute to 
a < 1. Experimental data indicate, that the transition from incoherent to coherent particle produc- 
tion happens at RHIC energies fovJ/y production. In the following we will discuss these trends in 
more detail and present calculations involving gluon shadowing and effects of energy-momentum 
conservation. 
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Figure 3: a vs. if (and xf) in proton-nucleus collisions at different energies for production of light particles 
(left) and J/y (right). Experimental data are taken from p^-p]]. 



2. Gribov inelastic shadowing 



In the relativistic Gribov theory [J|, |J] the collision proceeds through simultaneous interactions 
of the projectile with nucleons in the nucleus and therefore the intermediate states, shown in Fig. ||, 
are no longer the same as the initial state. The multiparticle content of the diagrams is given by 



AGK cutting rules [ 15 ], where the intermediate states are on-shell. The cut contribution of the dou- 
ble rescattering diagram can be expressed in terms of diffractive deep inelastic scattering (DDIS). 
The variable j8 = q^ m i = x/xp plays the same role for the Pomeron as the Bjorken variable, x, for 
the nucleon. We assume that the amplitude of the process is purely imaginary. The contribution 



from the second term in Eq. ( |1.2| ) to the total y*A cross section is given by 



y (2) 
'fA 



/i-M L 
d 2 b / dM 
J ML. 



danQ 2 ,x P ,r3) 



dM 2 dt 



\F A (q L ,b)\ 2 , (2.1) 



J t=o 



where 



dz p A {b,z)e iq L z 



(2.2) 



is the so-called longitudinal form factor and qi = ^—t m j n = m^xp 

In Eq. (2T), Mj nin corresponds to the minimal mass of the diffractively produced hadronic 
system, M^ nin = 4m 2 = 0.08 GeV 2 , and M} nax is chosen according to the condition: x P < xft"*. We 



use the standard choice for x'p ax 



0.1 Q17|]. It is convenient as it guarantees the disappearance of 
nuclear shadowing at x ~ 0.1 in accord to experimental data. 

We are interested in calculating shadowing for small-* quarks and gluons and consider be- 
cause of this DDIS on nucleons where this information can be extracted. Note that since Eq. ( p7T| ) 
is obtained under very general assumptions, i. e. analyticity and unitarity, it can be applied for 
arbitrary values of Q 2 provided x is very small 
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2.1 Diffractive gluon distribution 

The cross sections of diffractive processes are expressed through structure functions, which 
are in turn associated with distribution functions of partons in the Pomeron. The relation between 
the diffractive cross section and the diffractive structure function is given by 



do® N (Q 2 ,x P ,r3) 
dM 2 dt 



t=o 



4n 2 a em B P (3), n - P , 



Assuming Regge factorization we write the diffractive structure function as 



F^(x P ,Q 2 ,f3) = f P (xp)F(r3,Q 2 



(2.3) 



where the first factor is referred to as the (t -integrated) Pomeron flux and the second factor, F(j3 , Q 2 ), 
is the Pomeron structure function. 

Until recently, the diffractive structure function has been poorly known. In particular, the 
diffractive gluon density was affected by large uncertainty since it is not measured directly in ex- 
periment. Results of new high-precision measurements of the diffractive parton distribution func- 
tions (DPDFs) presented by the HI Collaboration [19, ^] give important constraints to our model 
of shadowing. Due to the indirect extraction of the diffractive gluon density, f}g J (f}, Q 2 ), from 
experimental data, two fits were presented, FIT A and FIT B, reflecting the systematic uncertainty 
of the procedure. Furthermore, a combined fit to DDIS data and diffractive di-jets [^Ij] results 
in a curve similar to FIT B yet with a slightly smaller gluon density. For completeness, we also 



compare the new results with the old HI parameterization Q22| ] presented in 2002. For details on 
extracted DPDFs and corresponding Pomeron parameters, C£>(0) and a P , we refer the reader to the 
experimental papers [ |T^ , 20]. 

The gluon distribution is almost a factor of 10 bigger than the quark distribution at the same 
Q 2 for a broad region of j8 fli9i 20]. Therefore, in the relevant kinematical range for hadron-nucleus 
and nucleus-nucleus collisions at RHIC and LHC, the gluon density dominates. In what follows, 



we will therefore consider structure functions of gluons in nuclei. The first term in Eq. ( |1.2| ) is 
then proportional to AGn{x,Q 2 ), while the second rescattering term is correspondingly equal to 

-A(A-\)G N (x,Q 2 ) Jd 2 bT 2 (b)f(x,Q 2 ), where 



f{x,Q 2 ) = 4% / dx P B(xp)fp(xp) 



Gn(x,Q 2 ) 



(fmin) 



(2.4) 



with B(xp) being the diffractive slope parameter. The gluon distribution of the nucleon, Gff(x, Q ) = 
xg(x, Q 2 ), was taken from the CTEQ6M parameterization [23]. The ratio under the integral in 
Eq. ( fO] ) is hence the ratio of gluon density in the Pomeron and in the nucleon. 



2.2 Models of multiple scattering 



The summation of all rescatterings in Eq. ( |1.2| ) is model dependent. In writing Eq. ( |2.4| ) we 
have assumed that the following factorization holds 



\F(q L ,b)\ 2 = Tl{b)Fl{t min ) 



(2.5) 
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Figure 4: Comparison of different models for higher order rescatterings. See text for details. 



where TaQj) = /_ ™ dzpA{b,z) is the nuclear density profile and Fa is given by 



F A (t, 



A \hnin I 



d 2 bJ (V^t~b)T A (b) 



(2.6) 



Here Jq(x) denotes the Bessel function of the first kind. Equation (2^5) is an identity for nuclear 
densities which depend separately on b and z, however we have checked that calculations with a 
Woods-Saxon nuclear density profile lead to negligible corrections to the exact expression. For 
consistency, we show calculations of (cr 1 ) + o^)/ Cy^ in Fig- % where for the dash-double-dotted 
curve Eq. ( |2~5| ) (denoted 'our' factorization) has been employed and the long-dashed curve shows 
the result using directly Eq. (2.1 ) (denoted 'KK' formula) [|16|]. The two curves practically coincide 
and indicate furthermore that higher order corrections are essential. 

We will now compute the total 7* A cross section using two models for higher-order rescatter- 



ings: a Schwimmer unitarization [ |24] ] which is obtained from a summation of fan-diagrams with 
triple-Pomeron interactions and a quasi-eikonal unitarization. Nuclear shadowing is studied in 
terms of the ratios of cross sections per nucleon for different nuclei, defined as 



R{A/B) = - 



B Cy*A 



A Oy*B ' 



(2.7) 



as a function of x, which can in turn be expressed via structure functions of the different nuclei. 
The simplest case is B=N, then 



R" g cn {A/N){x) 
Rf(A/N)(x) 



d 2 b 



d l b 



T A (b) 



\ + {A-\)f(x,Qi)T A (b) 



(2.8) 



2(A-l)f(x,Q 2 ] 



{1 -exp[-2(A-l)T A (b)f(x,Q 2 )]} (2.9) 



for the Schwimmer and quasi-eikonal models, respectively. In Eqs. (2.8) and (|2.9[), f(x, Q ) is given 



by Eq. (2.4). Both expressions (2.8) and (Ej), expanded to the first non-trivial order, reproduce 



the the second order rescattering result in Eq. (2.1). The gluon shadowing ratio R„ for x < 0.1 is 
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Figure 5: Gluon shadowing for Pb (lead) for different virtualities, Q 



shown in Fig. ^ where the calculations for Schwimmer (quasi-eikonal) model are depicted with a 
solid (dotted) curve. 

The quasi-eikonal model gives a stronger shadowing effect than the Schwimmer model. In our 
framework shadowing can also be studied as a function of the impact parameter b, given by 



R s g ch (A/N) (b) 



1 



Rf(A/N){b) 



l + (A-l)f(x,Q 2 )T A (b) ' 
1 



(2.10) 



2(A-l)T A (b)f(x,Q 2 



{1 -tx V [-2(A-l)T A (b)f{x,Q 2 )]} . (2.11) 



2.3 Results for gluon shadowing and validity of the model 

Summing up, in the framework of the Glauber-Gribov model the total y*A cross section can 
be calculated in a straightforward way provided the total y*N cross section and the differential 
cross section for diffractive production are known. In what follows, calculations are made with the 



Schwimmer unitarization in Eq. (2.8) employing both FIT A and FIT B for the gluon density of the 
Pomeron. 

The resulting gluon shadowing of the structure function of lead (Pb) is depicted in Fig. || for 
different virtualities, Q 2 . The shadowing is quite strong for x < 10~ 3 . We note that the QCD 
evolution of the main term and rescattering terms are effectively treated separately in our approach, 
and therefore the shadowing correction has a slow, logarithmic dependence on Q 2 [18]. Details and 
further results of the calculations can be found in [Eq]. 



Equation (2.8) does not take into account anti-shadowing effects which may play an important 
role for x>0. 1 . The model for diffractive production by virtual photon described above is applicable 
for intermediate Q 2 > 1 — 2 GeV 2 [|l^, ^(J as the parameterization of HI data leads to a violation 
of unitarity for low Q 2 and x — » and should be modified at very low x. This is clearly seen in 
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Figure 6: Pumplin ratio for different virtualities, Q 2 . The shaded area denotes the uncertainty of diffractive 
gluon distribution function from HERA JT^ , ^o| l. 



Fig. |H where we have plotted the so-called Pumplin ratio [25] 

jf x dx P B(x P )f P (x P )f3g^(f3,Q 2 ) 



R 



Gn(x,Q 2 



(2.12) 



which is in fact the ratio of diffractive and inclusive dijet production. Unitarity is violated when 
R> 1/2. Calculations show that this takes place at low Q 2 and low x. The shaded areas in Fig. || 
represent the uncertainty in the diffractive gluon distributions from FIT A and FIT B. We conclude 
that our calculations are reliable in the region x > 10~ 4 which is relevant for RHIC and most 
experiments at LHC, and should be taken with care for x < 10~ 4 at low Q 2 . 



3. Light hadron production at SPS and RHIC 

The model of gluon shadowing is now employed to study particle production in dAu collisions 
at RHIC energy = 200 GeV. There has been observed an increasing suppression of the nuclear 
modification factor (NMF) 



RdAu(p±,ri] 



1 d 2 N dAlt I djc^dT] 



<A^//} d 2 NZ/d P± dr] 



(3.1) 



with increasing pseudorapidity for charged particles [12]. We will not consider the effect of p±- 
broadening, or Cronin effect [p7|], at the moment, but rather study the r\ dependence of suppression. 
We assume, that in the ratio of forward to mid-rapidity nuclear modification factor 



R 



RdAu(p±,rj) 



(3.2) 
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Figure 7: Ratio of forward and mid-rapidity nuclear modification factor. Gluon shadowing contribution to 
the suppression in rj is depicted for two fits of diffractive gluon density. Data is taken from [jjJJ]. 



the Cronin effect is cancelled out (at this energy the effect is < 15%). In Eq. ( |3.2| ), R^Ai" * s tne 
nuclear modification factor at mid-rapidity divided by our calculations of gluon shadowing. The 
relation of kinematical variables is given by the standard formula 

x=^e-\ (3.3) 

where c ~ 3. In Fig. f7| we see that gluon shadowing contribute to the suppression of the nuclear 
modification factor at forward rapidities, but is not sufficiently strong to describe the data quanti- 



tatively. Equation (3.3) describes de facto mono-jet production. The mean Bjorken x involved in 



particle production calculated within perturbative QCD, which describes 2^2 particle collisions 



or so-called two-jet kinematics, is almost two orders of magnitude larger [28]. In this case, the 
effect of shadowing on the suppression in Fig. [7] would be strongly reduced. 

So far we have neglected a very important part of the nuclear suppression mechanism, namely 
the conservation of energy-momentum. This mechanism is responsible for the observed shape 
of a{xf) in the forward region as seen in Fig. [| Energy-momentum conservation affects the 



AGK cutting rules at finite energies []29J and results in an additional suppression factor []8J, [llj] 
corresponding to the substitution 

T A (b) — T A {b)tx V [-a e h l/{x +lPi _)AT A (b)} (3.4) 
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Figure 8: Comparison of nuclear suppression from pPb collisions at ^/s = 17.3 GeV (red points) [pOp, 
and dAu collisions at ^fs = 200 GeV (black points) [jOl [Hp, for two values of fixed xp . Red curve is 
extrapolation of low-energy suppression to high p±, green curve denotes the combined effect of shadowing 
and energy-momentum conservation at RHIC (solid: one-jet kinematics, dashed: two-jet kinematics). We 
present also predictions for pPb collisions at LHC energy (dash-dotted line). 

in the numerator of Eq. (j2"lf). The shape of o e JJ was formulated to describe data on low-p± nuclear 
modification factor for pPb collisions at SPS energy yfs = 17.3 GeV p0[], namely 



T eff 



(x+,p±) 



2 i 2 ' 

Pl+Po 



(3.5) 



and x + = 0.5 ( y x\ + Am\/s+xp) where m± is the transverse mass of the produced particle. Once 
the two free parameters, oq and p , are fitted to data at x F = 0, Eq. (p\4[) describe the SPS nuclear 



modification factor for < xp < 0.4. The same parameter values are also taken in the calculations 
of the suppression at RHIC. 

A comparison of SPS [ f30| ] (red circles) and RHIC [|l2|, 31] (black squares) data on nuclear 
suppression in pA collisions at xp = 0. 175 and xp = 0.375 together with calculations of the com- 
bined effect of energy-momentum conservation and gluon shadowing is presented in Fig. ||. The 
red dotted curve is an extrapolation of low-energy data to higher p± using Eq. (|3.4[). The green 



dashed curve show calculations using Eq. ( p.3| ) for one-jet kinematics. We have also calculated nu- 
clear effects at RHIC using the mean Bjorken x of a parton in the nucleus calculated within pQCD 
(solid green curve in Fig. |J). Predictions for LHC are given as well (dash-dotted curves). 

The most striking fact is that data at both energies seem to overlap. This indicates that a 
common mechanism dominate the suppression at forward rapidities at both top SPS and top RHIC 
energies. Moreover, since Gftf(x+ = 0) decrease with energy, RHIC data gives room for an ad- 
ditional small shadowing contribution. Unfortunately, our model give no information about the 
kinematics of particle production at low and moderate p±, and so the 'one-jet' and 'two-jet' curves 
indicate the total uncertainty of the model. At high-/?^, 'two-jet' kinematics is theoretically more 
justified. 



4. Nuclear effects in heavy quarkonium production 

Production of heavy state, such as Drell-Yan and heavy-flavor, give additional information on 
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the energy dependence of nuclear suppression mechanisms. Since leptons interact very weakly 
with the nuclear medium, Drell-Yan production holds information about the initial state effects 
in pA collisions, e. g. shadowing. On the other hand, hidden heavy-flavor is believed to interact 
quite strongly with the surrounding medium, either partonic or hadronic. Recently, the substantial 
decrease of the nuclear absorption in J / iff production in hA collisions between SPS energies, with 



Oabs ~ 4 mb []32|], and RHIC energies, with o a bs ~ 1 — 2 mb []13|], has attracted a lot of attention as 
it is in contradiction with the expectations of several theoretical models [§, [35|]. In what follows, 
we shall show that the apparent observation of the reduction of o a bs can be interpreted as a signal 
of the onset of coherent scattering for heavy state production. 

For heavy states the mass of the heavy system, M q q, introduces a new critical energy scale 

_ M 2 QQ R A m N 

SM-—^, (4.D 



It was shown in Ref. Qll| ] that AGK cutting rules are changed at s = sm- At energies below sm 
longitudinally ordered rescattering of the heavy system takes place and the produced heavy system 
is subject to nuclear absorption. In this energy interval Drell-Yan production does not experience 
nuclear suppression. At s > sm the heavy state in the projectile, which also includes light degrees 
of freedom, scatters coherently off nucleons of a nucleus and the conventional treatment of nuclear 
absorption is not adequate. For xp close to zero (central rapidity region) values of sm for J/yr 
and T belong to the RHIC energy region. In this kinematical region the effects of shadowing of 
nuclear partons become important. A similar approach for the description of // ^-suppression in 



dAu collisions at RHIC has been considered in Ref. [33]. For xp ~ 1 the high-energy regime 
corresponds to an interaction with a nucleus of the fluctuation of a projectile, containing heavy 
quarks (of the type of "intrinsic charm" mechanism of Refs. Jlo| , 34]). 



Thus at energies s < sm the inclusive cross section for production of particle a in a hA collision 
is given by 



E ^k {x+) _ w U '-tlW)^ ^) , (4 . 2) 

dp d> J t,{x+)a Q Q 



where R Sch is the ^-dependent shadowing factor calculated in Eq. ( |2.10[ ) and t, (x + ) = (1 — e) + ex 



8 



determines the x + dependence of absorption. The suppression is concentrated at much higher xp 
for QQ production than for the light hadrons because of the large mass of the QQ system, e. g. y = 2 



for J / y and y ~ 3 for T production. Equation (4.2) gives a good description of experimental data 



on charmonium production in pA collisions at £"lab^ 800 GeV/c with o q q = 20 mb and s = 0.75 



[11]. This corresponds to an absorption cross section of a a b s = 5 mb. Note that o q q is rather large, 
indicating that the cc pair is produced in the color octet state rather than in the colorless state. It 
can also be viewed as a DD (D*D*) system. 

Equation ( JO| ) is not valid at asymptotic energies as the assumption of longitudinal ordering 
leading to it is only valid at s < sm- For energies higher than sm the expression will change due to 
the correct treatment of coherence effects [|35|l 



1 - tx V [-B,{x + )o QQ AT A {b)} 



^{x + )o { 



QQ 



AT A (b)exp(-a e ^ f (xp)AT A (b)/2) , (4.3) 
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Figure 9: Centrality dependence of the nuclear modification factor in dAu (pPb) collisions at ^fs = 200 
GeV (5.5 TeV) for J/y at different rapidities. Data are taken from j 



which is similar to the energy-momentum conservation effect for light quarks. In the model of 
Ref. [35], Gqq was found to be the QQ — N total cross section. 

We would like to argue that this leads to an incorrect behaviour at high energies, and propose 
an alternative procedure. Considering non-enhanced, Glauber-type diagrams the effective cross 



section cfJJ should be proportional to x r + , thus satisfying AGK cancellation in this limit. It was 



shown in Ref. Ql 1| ] that at xj? ~ 1 the second rescattering in the low and high energy limits should 
coincide. This means that a^J £x\o Q q. Experiment on J/y production in dAu collisions at 
RHIC [13] was performed in the central rapidity region, where x + ~ 0.025 — 0.05 and o e Q ^ is 
therefore very small. This favors the pure nuclear shadowing scenario, which means that we do not 
include contributions from non-enhanced, Glauber-type diagrams. 

Results of our calculations for J j\\f and T production in dAu collisions at RHIC are shown 
in Figs. H and 10 together with data from [13] (solid curves). Both the centrality and the rapidity 
dependence of the data is well reproduced. We predict a stronger shadowing effect in pPb collisions 
at LHC energy ^fs = 5.5 TeV, which is given by the dashed curves in Figs. ^ and 1C. This effect 
should also be taken into account when modeling nucleus-nucleus collisions at high energy. 



5. Conclusions 

We have calculated gluon shadowing in the Glauber-Gribov model using latest parameteriza- 
tion of diffractive gluon distribution from HERA. A strong shadowing effect is found. Our model is 
applicable for Q 2 > 2 GeV 2 and 10 5 < x < 0.1, well suited for analysis of moderate and high-/?^ 
particle as well as heavy-flavor production at high-energy experiments. 
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Figure 10: Rapidity dependence of the nuclear modification factor in min. bias dAu collisions (pPb) at 



*fs = 200 GeV (5.5 TeV) for J / y and T at different rapidities. Data are taken from [ 13 1. 



Particle production at RHIC is dominated by coherent production both for light and heavy 
particles at low and moderate p±. This is most clearly observed at mid-rapidity where a(xp = 0) 
is below but close to unity for both charged particles and J/y - the small suppression is solely due 
to gluon shadowing. Shadowing effects are stronger for light than for heavy particles as expected. 

Both centrality and rapidity dependence of J/y production in dAu collisions at RHIC have 
been described with our model. Even so, gluon shadowing alone cannot explain the r\ dependent 
suppression of light particles measured at RHIC. At forward rapidities energy-momentum conser- 
vation comes into play and contributes strongly to the observed suppression. This happens more 
rapidly for light particle production, for J/\j/ and T this effect is shifted to larger values of xp due 
to the large mass. The combined effect of energy-momentum conservation and gluon shadowing 
shows good agreement with SPS and RHIC data. 

The discussion of mid-rapidity nuclear modification factor is out of scope of this paper. A 
detailed calculation with the above mentioned effects and Cronin effect should be performed at 
all rapidities to draw quantitative conclusions from experimental data. Gluon shadowing by itself 
can also be checked against dilepton or direct photon data, and also in ultra-peripheral heavy-ion 
collisions. 
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